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2 : ABSTRACT 

We present optical and near-IR (OIR) observations of the major outbursts 
of the neutron star soft X-ray transient binary system Aquila X-1, from summer 
K^ \ 1998 - fall 2007. The major outbursts of the source over the observed timespan 

T^-j- \ seem to exhibit two main types of light curve morphologies, (a) the classical 

rn 



Fast-Rise and Exponential- Decay (FRED) type outburst seen in many soft X-ray 

transients and (b) the Low-Intensity State (LIS) where the optical-to-soft-X-ray 

O ' flux ratio is much higher than that seen during a FRED. Thus there is no single 

Q \ correlation between the optical (R-band) and soft X-ray (1.5-12 keV, as seen 

by the ASM onboard RXTE) fluxes even within the hard state for Aquila X- 

\rl ' 1, suggesting that LISs and FREDs have fundamentally different accretion flow 

Vh ' properties. Time evolution of the OIR fluxes during the major LIS and FRED 

outbursts is compatible with thermal heating of the irradiated outer accretion 

disk. No signature of X-ray spectral state changes or any compact jet are seen 

in the OIR, showing that the OIR color-magnitude diagram (CMD) can be used 

as a diagnostic tool to separate thermal and non-thermal radiation from X-ray 

binaries where orbital and physical parameters of the system are reasonably well 

known. We suggest that the LIS may be caused by truncation of the inner disk 

in a relatively high M state, possibly due to matter being diverted into a weak 

outflow. 

Subject headings: accretion, accretion disks — stars: neutron — stars: — X-rays: 
binaries — individual (Aquila X-1) 
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Introduction 



Aquila X-1 is one of the most frequently recurring soft X-ray transients in the sky. It 
goes into outburst approximately once every year although the periodicity is not fixed (ISimon 
20021 ). The frequent outburst recurrence makes it an ideal source for studying accretion 
flow during an outburst as well as changes in accretion flow from outburst to outburst. 
The physical and orbital parameters of this system are relatively well known . The orbit a l 
period of the binary system was determined to be 18.71 ± 0. 6 hou rs by IWelsh et al.l (J2000l ) 
from optical photometry during quiescence. IChevalier et al.l (119991 ) estimated the quiescent 
V-band magnitude to be 21.6 ± 0.10 magnitude. However, there is a con taminating star 



abou t 0.48" East of Aquila X-1 which has a V magnitude of 19.42 ± 0.06 (IChevalier et al. 

19991 ). The distance to the system is between 4-6.5 kpc as inferred both from (a) limits 

)e flll ing donor star of spectral type between K6-M0 
20011) and fb) pea k X-ray fluxes during type-1 bursts 



on optical brightness of the Roche lo 



Chevalier et al 



1999 



Rutledge et al 



( ICzerny. Czerny. fc Grindlay 



198?! : iRutledge et al.ll200ll ) 



In systems like Aquila X-1 and other Roche lobe overflowing so ft X-ray trans ient (SXT) 
binary systems, the thermal viscous disk instability model (DIM; iLasotal 120011 ) is usually 
invoked to explain the typical fast-ri se and exponential- decay (FRED) type outburst light 
curve morphology (jChen et al.l 119971 ). The DIM alone however cannot produce realistic 
FRED outbursts. Another mechanism that plays a key role during outbursts of SXTs is 
X-ray irradiation. The irradiation is caused by the intense X-rays from the inner regions of 
the accretion disk and neutron star surface (if the system harbors a neutron star) imping- 



i ng on the outer di s 



( Cunningham 



1995 



1976 



Dubus et al. 



<:/donor, and by r e -emission of this impinged 



Vrtilek et al. 



199 



i 



1990 



O'Brien et al 



energy 



van Paradiis fc McCliiitock 



20021 : iHvnes et al 



2002 



1994 



in lower energies 



McChntock et al. 



Russell et al.ll2006[ ). Another 



catagory of modiflcations to the DIM is generated by the possible truncation of the inner 
region of the accretion flow, where the disk inflow can be replaced by a radiativ ely ineffi- 
cient accretion flow (RIAF). There is an increasing body of work on th i s topic (llchimaru 
I977I : iRees et al.lll982l : [Naravan fc Yilll994l : Isiandford fc Begelmanlll999l : iMeieilboOlh . but 
the overall impact of these efforts on the outburst morphology is currently still unclear. 

CoUimated relativisti c outflows, lo osely termed as "jets", have been observed in several 
soft X-ray transients (see iFenderl |2006| . for a review) and their broad band spectral energy 
distribution (SED) from radio to X-rays have been fltted with jet-fdisk models comprising 

of a freely expanding "Bla ndford - Konigl" jet (JBlandford fc Konigllll979l : lFalcke fc BiermannI 

1999 : lMarkoff et al.ll2005l ) and a thermal- viscous "Shakura-Sunyaev" accretion disk ( Shakura fc Sunyaev 
I973I ). Observations suggest that compared to black hole (BH) systems, the neutron star 
(NS) SXTs tend to have intrinsically weaker jets and also the transition from optically-thick 
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to optically-thin jet spect rum in NS systems probably occurs in the mid-IR (iMigliari &: Fender 



2002 



2006 : Mighari et al 



2006 



this work), compared to that in near-IRfor B H systems (I Corbel fc Fender 



Buxton fc Bailvnll2004J : iHoman et al.ll2005l : iRussell et all 120061 ) 



During outburst, light curve morphology of many sources are not FREDs. The X-ray 
spectral and temporal features also show an enormous di versity during outbursts of SXT 
syst ems, and are categjorize d into various states. See e.g. iMcClintock fc RemillardI (20061) 
andlHoman fc Bellonil (120051 ) for extensive discussions on definition of states and Ivan der Klis 
( 20061 ) for their timing properties. Despite the multitude of states observed in the different 
systems and different outbursts, two "canonical" states can be identified during the outburst 
of most SXTs, characterised by presence or absence of a thermal and nonthermal component 
in the energy spectrum. In the temporal domain, presence or absence of a thermal component 
usually shows up respectively as low or high rms variability in the light curve. A host of other 
intermediate and extreme states are also seen, which could be a superposition of the thermal 
and nonthermal components with varying contributions. Precise determination of X-ray 
state requires modeling of spectral energy distribution and Fourier power density spectra. 
Since this work primarily focusses mainly on OIR observations, we have used dates of X- 
ray state transitions of A quila X-1 from published st ud ies of X-ray spec t ral and temporal 
prop e rties of the source (Maccarone &: Coppil l2003a| ]bl: iReig et al.l 12004 : iMaitra fc Bailyn 
2004 : [Rodriguez et al.l l2006l : IYu fc Dolencd 120071 ) to study correlation between X-ray state 
and OIR fiux-density. 

We have observed Aquila X-1 for ~ 9 years, through many outbursts, to explore the 
nature of outburst morphologies in a single object. Even though the OIR (as well as X-ray) 
light curve morphologies from outburst to outburst differ, we find systematic and repeated 
tracks traced by the source in OIR color-magnitude diagrams. Since we have long-term data 
in two bandpasses (R and J) only, we have used a simple model to show that the color 
evolution and the resulting track are consistent with heating of a constant area blackbody 
and can be caused by an irradiated outer disc. 

In [|2] we describe the general reduction procedures adopted in this work. The light 
curve morphologies in J, R and ASM bands (of the RXTE satellite) and optical-soft X-ray 
correlations during the major outbursts are described in §31 Evolution of Aquila X-1 on the 
OIR color-magnitude diagram is presented in ^ Finally we present our discussions in ^ 



-4- 

2. Observations and Data Analysis 

2.1. Optical and Near-IR data 

Our long term monitoring campaign used the Yale 1.0m telescope at Cerro Tololo Inter- 
American Observatory (CTIO) in Chile to observe Aquila X-1 from 1998-2002 on a daily 
basis (weather and Sun-angle constraints permitting). Since February 2003, the CTIO 1.3m 
telescope has been used for the monitoring. In both cases, images were obtained using the 
ANDICAMj instrument, which is a dual-channel imager capable of simultaneously recording 
images at an optical and an infrared wavelength. Prior to MJD 52430 (2002 June 5), the 
optical R-band data were recorded by a Lick/Loral-3 2048 x 2048 CCD. Thereafter, the 
R-band has been recorded by a Fairchild 447 2048 x 2048 CCD. The near-IR J-band data 
are recorded by a Rockwell 1024 x 1024 HgCdTe Astronomical Wide Area Infrared Imager 
(HAWAII) imager. Between MJD 51130-51367 (1998 November 13 - 1999 July 8th), the 
source was observed using a "wide-R" filter. The wide-R filter has a broader bandpass than 
the Johnson R filter used during the rest of the campaign. Even though a determination 
of the Johnson R-band flux cannot be made from the wide-R-band measurements, the light 
curve morphology changes little due to the change of filter. We calculated the quiescent 
wide-R-band differential magnitude (between MJD 51232 and MJD 51306) and applied an 
offset to the differential magnitudes obtained using wide-R filter so that the quiescent level 
was set to 18.8 magnitude, which is the quiescent R-band magnitude of Aquila X-1. 

Under normal viewing conditions at CTIO, it is extremely difficult to separate the 
true optical counterpart of Aquila X-1 from the bright contaminating star. Therefore the 
magnitudes reported in this work are the combined magnitudes of the Aquila X-1 system 
and the contaminating star, and when we say magnitude of Aquila X-1 we actually mean the 
combined magnitude of Aquila X-1 and the contaminating star. In quiescence, light from 
the contaminating star dominates, but once an outburst begins, photons from Aquila X-1 
dominate. 

The optical data reduction, including bias subtraction, overscan correction, and flat- 
fielding, was done through the standard IRAF data reduction pipeline developed for these 
instrumental. For the J-band near-IR data, multiple dithered frames were taken and then 
fiat-fielded, sky subtracted, aligned and average-combined using an in- house IRAF script. 
During a large portion of the year 2004, one of the four quadrants of the IR CCD was 
not working. In most of the observations during this period, one or more of the reference 



^http://www. astro. yale.edu/smarts/ANDIC AM 

^For details see |http://www. astro. yale.edu/smarts/smartsl3m/optprocessing.txt 
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stars were not present in the frames. We used whichever comparison stars were present and 
adjusted the resulting differential magnitudes appropriately. 

We used the field stars 2MASS 19111739+0035039 and 2MASS 19111633+0034333 with 
J-band magnitudes of 12.404 and 12.490 respectiveljo, to obtain differential magnitude of the 
source. Flux calibration of these reference stars in the R-band was performed by comparing 
with Landolt standard stars on the photometric night of MJD 53904 (2006 June 18). In 
order to estimate time variability in flux from any of the reference stars, we calculated the 
difference between the observed magnitudes of the reference stars on every night. Smaller 
than 0.02 mag in R-band, this difference is the primary source of error in flux determination 
because during outbursts the Poisson error due to photon count-rate is much smaller. The 
S/N is even better in J-band where interstellar extinction is lesser. 



For magnitude to flux conversion, following iBessell et al.l (119981 ) and ICampins et al. 



(1l985l ). we assumed that R=0 magnitude converts to 3064 Jy and J=0 magnitude converts 
to 1603 Jy. IWelsh et al.l (120001 ) estimated that 83 ± 2% of the quiescent R-band flux is from 
the contaminating star. Assuming that the brightness of the contaminating star is constant 
in time, we subtracted this flux from the total observed flux to obtain the flux of Aquila X- 
1 only, during quiescence as well as outbursts. Interstellar reddening and extinction was 



calculated ass uming a hydrogen eg 



of Aquila X-1 (JDickey fc Lockman 



umn density [Nh) of 3.4 x 10^^ atoms/cm^ in the direction 



19901). From deep optical photometry and spectroscopy of 



the source during quiescence IChevalier et al.l (119991) estimated the color excess E{B — V) = 
0.5 + 0.1 which is consistent with the lDickey fc LockmanI (Il990l ) value of Nh in the direction 
of Aquila X-1. For dereddening the OIR flux we used the standar d A^^ to visual extinction 
{Ay) correlation given hy Nm/ Ay = 1.79 x 10^ ^ atoms/cm^/mag (jPredehl fc Schmittlll995l ) 
and the extinction law of ICardelli et al.l (119891 ) which gives Ar/Av = 0.748 magnitude and 
Aj/Ay = 0.282 magnitude for Aquila X-1. Taking reddening and extinction in the direction 
of Aquila X-1 into account, the observed R — J color is linearly related to the corresponding 
intrinsic spectral index (aju-int) by the relation 

logifn) - log{f/ 



O^RJ-Ant 



logiuRj - log{uj) 



2.47 - 1.63 x(i?-J)o,„ 



erved- 



2.2. X-ray data 



The X-ray data was taken from the A ll Sky Monitor (AS M; iLevine et al.lll996l ) aboard 
the Rossi X-ray Timing Explorer (RXTE; iJahoda et al.lll996f ) satellite. The ASM has been 



^Taken from the 2MASS point source catalogue at lhttp://www. ipac.caltech.edu/2mass 
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monitoring the sky in the 1.5-3, 3-5 and 5-12 keV energy bands continuously since 1996. 
Usually a given source is observed a few times every day, subject to constraints like angular 
proximity to the Sun etc. The results are made publicly available by the MIT/ASM team 
through their websito We have used the one day average light curve for Aquila X-1. During 
quiescence, the flux from Aquila X-1 is below the detection limit of the ASM, however since 
we are interested only in the major outbursts, which are detectable by the ASM, we did 
not take into account data during quiescence, for our optical-X-ray correlation study. Also, 
during the early as well as the late stages of the outbursts, the count-rates in the individual 
energy channels are small; therefore we have only looked at the summed count-rates in this 
work. The ASM count-rates were converted to flux-density at 5.2 keV, assuming a Crab-like 
spectrum. The Crab spectrum in the ASM energy range is a simple power law with a photon - 



index of -2.1 and the total flux within that range is 2.8 x 10~^ ergs/cm^/s (jSmith et al.ll2002l ). 
so the fiux-density at 5.2 keV is 1.06 mJy. The flux conversion (from ASM count-rate to /^Jy) 
is therefore a numerical factor by which the total ASM count rate is multiplied. The choice 
of 5.2 keV was motivated by the fact that it is the geometric mean of the half-transmission 
points of the ASM response shapqj near 2.3 and 12 keV. We used geometric mean (CM) 
because the GM compensates for the fact that almost all X-ray spectra slope downward in 
count rate toward higher energies. 

Evolution of spectral shape during an outburst introduces uncertainities in flux de- 
termination. Pivoting uncertainities (due to spectral evolution) in ASM's X-ray range are 
small compared to the optical to X-ray separation in frequency space. We used the online 
WebPIMMS^ tool which shows that for a given total ASM count-rate, the difference in flux 
in the ASM band-pass is less than 10% between a power law of photon index 2.1 (i.e. Crab 
spectrum), a power law of photon index 1.7 or a blackbody at 2 ke V (typical photon in 



dex/blackbody temperature of Aquila X-1 in hard / soft state; see e.g. iMaccarone fc Coppi 



2003bl : iMaitra fc Bailynll2004j : [Rodriguez et al.ll2006l ). Therefore we added a 10% systematic 



error, in quadrature to the Poisson error, during the ASM count-rate to flux conversion to 
account for the uncertainity in the spectral shape. The precise value for pinning and rescal- 
ing the flux conversion is not very sensitive to the choice of keV reference, as long as this is 
done self-consistently (i.e. scale to Crab flux density at the reference keV, and implemented 
in this work). 



* http://xtc.mit.edu 



^Taken from appendix F of the RXTE technical appendix (http://heasarc.nasa.gov/docs/xte/appendix_f.htmll 
and Ron Rcmillard, priv. comm. 

^http://heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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In this process, any information about changes in the energy spectrum is lost, but we are 
primarily concerned with overall X-ray flux here, and using all the counts results in better 
signal-to-noise ratio. In order to assess the X-ray state of the source, we ha ye used the d ates o f 



state transitions f rom p u blished literature iMaccarone fc Coppil ( 2003a Ibf) ; iReig et al.l (120041 ) ; 
Maitra fc Bailvnl J2004J ): [Rodriguez eTaP J2006h : lYu fc Dolencell2007h . For the most recent 



outbursts (e.g. outburst #9, #10) with no published dates of state transition(s), we used 
the (9.7-16 keV/6-9.7 keV) hardness ratio from archival RXTE/PCA pointed observations 
(Manuel Linares, priv. comm.) whose value correlates strongly with the X-ray state of the 
source. For outbursts with sparse coverage of pointed observations (e.g. outburst #5) we 
used the (5-12 keV/3-5 keV) hardness ratio from the RXTE/ASM data to infer the spectral 
state. 



3. Major Outbursts of Aquila X-1 since 1998 June 20 

In Fig. [T]we show the results of long-term multi-wavelength monitoring of Aquila X-1. 
The RXTE/ASM monitoring of Aquila X-1 dates back to the beginning of 1996, but our 
nightly R-band monitoring started on MJD 50984 (1998 June 20th). The J-band observations 
started MJD 51677 (2000 May 13th) onwards, but the coverage was somewhat sporadic till 
the end of 2002, as can be seen from the top panel of Fig. [H During quiescent periods, the 
system is detected in both of the OIR bands but not in ASM. Therefore in Fig. [T] we have 
only plotted those ASM points for which the signal-to-noise ratio is better than four. 

Since observing constraints sometimes make it difficult to date the precise starting and 
ending moments of an outburst, in Tabled] we note the MJD when the ASM, R-band and 
J-band fluxes attained a maximum during an outburst. We define a major outburst as a 
period of activity of the source, spanning at least 20 days or longer, when the source fiux 
was at least 3a above quiescence in at least one of the observed (J, R and ASM) bandpasses. 
Eleven major periods of activity were observed by the ASM during this period. All these 
periods of major activity were also seen in OIR, except the one during late 2005-early 2006 
when the source's angular proximity to the Sun was too close to observe with ground based 
OIR telescopes. Even the RXTE/ASM data, while showing unmistakeable signs of activity 
during this particular outburst, does not have sufficient coverage and S/N to classify the 
light curve morphology. 

Besides the major outbursts, several smaller and shorter periods of activity are seen in 
the OIR light curve (Fig. [1]). Because of the low signal-to- noise ratio, particularly in ASM 
data, and scarcity of PGA pointings during these periods of minor activity, they have not 
been discussed in the present work. However, the time of peak OIR flux from the source 



during these "mini" outbursts have been presented in Table [H 



3.1. Lightcurve morphology: FREDs and LISs 

From the long-term OIR monitoring of Aquila X-1 from MJD 50984 - MJD 54411 (1998 
June 20 - 2007 Nov 7) only three out of the major OIR outbursts are purely FREDs. The 
remaining are a combination of FREDs with another active state where the flux is above 
quiescence, but does not rise or decay on the timescale of a FRED. Instead, a plateau type 
light curve with significant variability on a timescale of days is seen in all bands. Sometimes 
this state persists significantly longer than the FREDs. For Aquila X-1 the typical observed 
FRED rise times are ~ few days to a week and the decay times are ~ 3-4 weeks. On the 
other han d, the longer l a sting, somewhat variable state can last for longer than a month. 



Following IWachter et al.l (120021 ). who reported observation of such state in the neutron star 
system X1608-52, we call this a "Low-Intensity State" (LIS). In Fig. [2]- Fig. Owe show 
blowups of the R-band and ASM light curves for each of the major outbursts in Fig. [T] 
(also tabulated in Table [1]). FREDs and LISs are shown by different symbols in Fig. [2]- 
Fig. [5l During a LIS, the OIR flux is near the maximum, but always slightly lower than 
the corresponding maximum flux observed during a pure FRED outburst. For LISs, the 
variability in the light curve is more complicated than the typical 'up to a maximum, then 
down to quiescence' morphology seen for FREDs. For a comparable optical brightness, the 
soft X-ray flux is a factor of seven or more lower in a LIS than in a FRED. This high optical- 
to- X-ray flux-ratio, discussed in greater detail in §3.2[ is one of the defining characteristics 
of a LIS. 



3.2. Optical— soft X-ray correlations 

The strongest evidence for the existence of an independent LIS comes from optical-X-ray 
correlations. From the long-term OIR monitoring data, we estimated R-band fiux-density 
from Aquila X-1 as described in §2.1[ The RXTE/ASM one day averaged count-rates were 
converted to fiux-density as well ( §2.2p . Both light curves were then searched for quasi- 
simultaneous optical-X-ray data where the observation times in each bandpass were not 
separated by more than one day. In Fig. [6] we show the optical-soft-X-ray fiux correlation 
during the major outbursts for which we had quasi-simultaneous ASM and R-band data. 

The optical and soft-X-rays increase approximately linearly with f opt/ f x-ray ~ 1.3 — 1.5 
for FRED outbursts and ~ 9 — 12 for the LISs. Given the considerable scatter in the data. 
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we have not attempted to fit the data with any function and the fiux-ratios quoted above, as 
shown in Fig. [6|, are a guide to the eye. For the same X-ray brightness, there is no strong hint 
of hysteresis, i.e. the optical brightness during outburst rise is not significantly different from 
than that during the outburst decline. The absence of optical-X-ray hysteresis in Aquila X-1, 
goes to support the hypothesis that either jets in neutron star systems are generally weaker 



or the jet break oc curs at a lower frequency (JMigliari fc Fended l2006l : iMigliari et al.l 12006 



Russell et all 120071 ). 



Both FREDs and LISs start/end the outburst at quiescence and hence have a common 
start /end-point near the lower-left region of fopt — f x-ray plane. For some of the weakest 
outbursts, the ASM fluxes are within a few sigma of detection threshold. At such low 
luminosities, both FREDs and LISs occupy similar regions on the fopt — f x-ray plane. It 
has been observed that for some low luminosity outbursts, the sour ce may never come ou t 



of the X-ray hard state during the entire period of activity (e.g. see [Rodriguez et al.ll2006l ). 
Even though the peak optical brightnesses of the brightest LISs are comparable to the peak 
optical brightness during FREDs, LISs are typically fainter than FREDs. The distinction 
between LIS and FRED outbursts stands out clearly when optical and X-ray data from all 
the major outbursts since 1999 June 20 are plotted on fopt — f x-ray plane (Fig. E]). LISs and 
FREDs occupy separate regions of the fopt — f x-ray plane and evolve on different timescales, 
suggesting different physical mechanisms for their origin. 



4. OIR color evolution 

In Fig. Hand Fig. [H] the time evolution of the R and J-band magnitudes and that of the 
R— J color is shown for the major outbursts that were observed in both bandpasses. There 
seems to be no significant inter-band time lag during rise or decay of an outburst. Irrespective 
of the light curve morphology during the outburst, the R— J color is always bluer during an 
outburst compared to the quiescent value, and therefore {R — J)outburst < {R — J)quiescence 
as seen in Fig. [7] and Fig. [HI This decrease in R— J color (becoming bluer) is characteristic 
of increased irradiation heating of the outer disk. 

During quiescence, the disk flux is very low compared to the combined flux from the 
contaminating star and the donor star. It is however interesting to note that the R— J color 
during an LIS is significantly redder than that during the peak of a FRED outburst. 

Fig. [9] and Fig. [10] shows the OIR color-magnitude diagram for the different outbursts. 
The R-band magnitude is plotted along the ordinate. As is typical for conventional color- 
magnitude diagrams, the ordinate axis has been reversed so that higher luminosities (i.e. 
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lower magnitude) are plotted near the top and vice- versa. The R— J color is plotted along 
the abscissa. 

During periods of quiescence, the source flux is low and essentially dominated by photons 
from the contaminating star and the donor, and it lingers near the bottom-right corner of 
the color-magnitude diagram. As the source enters an outburst, the brightness increases 
and the color becomes bluer, due to an increasing contribution of the accretion disk flux. In 
Fig. M and Fig. (TDl the observations when the source was in the hard state during the rise of 
an outburst, when it was in the hard state during the decay and when it was in the soft state 
are shown by different symbols. As noted in § 13.21 there a re also outbursts when the source 



spent the entire period of activity in an X-ray hard state (iRodriguez et al.ll2006l ). The OIR 
data during such outbursts are shown by open squares. There seems to be no significant 
change in OIR color across the X-ray spectral state transitions. 

We attempted to test the simple ansatz that the observed evolution on the color- 
magnitude diagram is consistent with the thermal heating of a blackbody. Our goal is 
to see whether the flux and color variations are correlated as would be expected if the only 
change is in the disk temperature. For simplicity, and because we had only two band-passes, 
we considered the entire annulus to be a single temperature blackbody, even though real 
accretion disks are likely to have a significant temperature gradient. In this simple model, 
the R— J color is determined by the temperature of the annulus and the R-band brightness 
is proportional to the projected area of the annulus and inversely proportional to the square 
of the distance to the source. Since the observed brightness depends on the projected area, 
we multiplied the area of a disk of radius Rdisk, by a multiplicative factor / < 1, and used 
this area as the area of the annulus. Thus in this simple model the brightness is propor- 
tional to f {Rdisk / D diskY cos (i) , where Ddisk and i are the source distance and i is the orbital 
inclination respectively. The only adjustable parameter here is /, which changes the overall 
projected area of the annulus. The parameter / is changed from outburst to outburst, keep- 
ing other p arameters constant. F itting the ellipsoidal variations observed in the quiescent 



light curve, iRobinson et al.l (120011 ) estimated that 36° < i < 55°. We assumed an inclination 
of 45°. An upper limit on Rdisk can be presumed to be equal to the radius of the Roche 
l obe of the accr etor and therefore can be estimated from the mass ratio and orbital period 



( lEggletonlll983l ). Evolution of the source on the color- magnitude diagram was generated by 
calculating the color and brightness of the annulus at varying temperature. We note that 
the track generated by the heating/cooling of this constant area blackbody broadly describes 
the evolution of the source during its outbursts. For a purely viscous disk (no irradiation) 
the spectral index between R and J bands (as defined in Eq. [T]) should lie between 1/3 and 
2. The observed range in i? — J color during outbursts is approximately between 1.35 and 
2.1, which using Eq. [1] translates to a spectral slope between 0.27 and —0.95 respectively. 
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The observed range in spectral indices is quite different from what could be expected from 
a purely viscous disk and favors irradiation heating. 



We found that even for the largest possible disk radius of ~5 light seconds (M, 



donor 



O.SMfT,, M 



05 



'■accretor 



1.5Mq, Porb = 18.95 hour) , the model-predicted R-band flux is too small 



for the lowest distance estimate of 4 kpc by iRutledge et al.l (120011 ). Using Ddisk = 3 kpc 
seems to lie close to the observed datapoints (see Fig. [9] and Fig. [10]). Since we are not doing 
detailed spectral fits, we do not claim the distance to the source to be 3 kpc, but show that 
this very simple model can roughly predict the outburst behavior of the accretion disk. The 
model curves overplotted in Fig. [9] and Fig. [10] are not fits to the datapoints but created 
using a disk of radius 5 light-seconds, inclined at an angle of 45° to the observer and at a 
distance of 3 kpc from the observer, viz. a set of parameters that is generally consistent with 
the observed values, and with a small range (0.6-0.8) in the value of the multiplicative factor 
/. The general consistency of the data over many outbursts supports the idea that the OIR 
flux is largely due to thermal emission from an irradiated disk. This is in contrast to recent 
observations of accreting black h ole systems in which a strongly non-thermal component 



is oft en observed in the near IR (jBuxton fc Bailyn] l2004l : JHoman et al.l l2005l : JRussell et al. 



20061 . Russell, Maitra & Fender, in prep.). 



5. Discussion 



Multi-wavelength monitoring of the neutron star soft X-ray transient Aquila X-1 in 
near-IR, optical and soft X-ray data from RXTE/ASM over the past nine years has allowed 
us to study the evolution of the source with unprecedented coverage. The OIR and X- 
ray light curves show a wide range of outburst light curve patterns. However the major 
outburst s show two distin ct light curve morphologies. The fast-rise and exponential- decay 
(FRED; IChen et al.lll997] ) type outbursts characteristic of many soft X-ray transients is seen 
in Aquila X-1 too. However, only three of the major outbursts seen since 1998 June show just 
a FRED ty pe light curve profile . More often a prolonged, somewhat variable "Low- Intensity 
state" (LIS: IWachter et al.ll2002l ). uniquely characterized by an extremely high optical-to-soft 
X-ray fiux-ratio, is seen either before or after a FRED. It has been observed that for some of 
these "soft X-ray faint" LI S outbursts, Aquila X -1 never comes out of the X-ray hard state 
during the entire outburst (JRodriguez et al.ll2006l ). Other than these major outbursts, small- 
scale flaring activities were also observed. While the origin of such small-scale short activity 
remains unknown, such a ctivity has been observed in other sources like XTE J1550-564 by 
Sturner &: Shraderl (J2005[) who suggested that these small outbursts could be due to discrete 
accretion events. We suggest the following framework for understanding these results: 
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Irradiation dominates OIR flux for hath FRED and LIS outbursts of Aquila X-1. During 
the course of FRED and LIS outbursts of Aquila X-1 the evolution of the source on the OIR 
color-magnitude diagram is consistent with heating/cooling of an irradiated outer accretion 
disk as a constant area blackbody. A classical viscous accretion disk (e.g. a Shakura- 
Sunyaev disk) spectrum has a frequency dependence of Fy ~ u^/^ between kT(Rout)/h <^ 
z/ <C kT{Rin)/h where T{Rin) and T{Rout) are the disc temperatures at the inner and outer 
edges of the disk. At frequencies much smaller than kT{Rout)/h the disk spectrum is simply 
the Rayleigh- Jeans region of a blackbody with F^, ~ z/^. For soft X-ray transients kT{Rin)/h 
is typically in far-UV or soft X-rays. Therefore in our case the flux ratios (a.k.a. colors) 
of two bands, from a purely viscous disc are expected to be constant as long as both bands 
are either on the viscous (u^^^) or Rayleigh- Jeans (z/^) regime. In the special case when 
the break frequency between the viscous dissipation and the Rayleigh- Jeans regime given 
by Inbreak = kT{Rout)/h lics between R and J bands, the spectral index should lie between 
1/3 and 2. The observed systematic and continuous change in color (and the corresponding 
range in the spectral index from —0.95 to 0.27; Q seen in Fig. [7] and Fig. [8] is strongly 
suggestive of an irradiative origin for the observed OIR emission and not of a purely viscous 
origin. 

Neither the OIR color nor its brightness change sharply during an X-ray spectral state 
transition. This confirms that for Aquila X-1 the outer accretion disk, which emits most of 
the thermal OIR radiation, is not affected by X-ray spectral state transitions which change 
the physical environment of the inner regions of the a ccretion disk. This is in s harp contrast 



with black hole binaries like 4U 1543-47, GX 339-4 flBuxton fc Bailvnl l2004l : iHoman et al. 



20051 ) where near IR colors are strongly correlated with X-ray state changes and the cor- 
relation interpreted as evidence for non-thermal emission. In FREDs, a slow increase in 
brightness above that predicted by our constant area blackbody model, as the source be- 
comes bluer and brighter, is most likely due to increase in the surface area of the disk (e.g 
most prominent for outburst #2 and #5 in Fig. [H]). It is possible that at highest mass accre- 
tion rates the high X-ray flux from the central inner disk causes the flaring of the outer disk 
and hence increase the effective surface area. In contrast, during the LISs, when the source 
probably spends the whole period of activity in hard X-ray state, there does not seem to be 
any signiflcant change in the disk area (e.g. outburst #7 in Fig. [TU]). Therefore it is likely 
that the effect of X-ray state transition on the outer accretion disk is both an increase of 
the outer disk temperature due to enhanced irradiation from the inner disk and an increase 
in the emitting area of the outer disk due to flaring. Total irradiation (hence OIR flux) is a 
bit smaller in LIS than in peak FRED. This could be due either to lower M, or to radiative 
inefficiency lowering L/M ratio. 

Broad band quasi- simultaneous observations in LIS suggest outflow. Soft X-ray flux is 
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very low in LIS, but hard X-rays are high (JMolkov et aLll2004j . Linares et al. in preparation,). 
While absence of soft, thermal X-ray photons suggests that the inner edge of the accretion 
disk is cold and probably receded during the LIS, the high OIR brightness (see Fig. for the 
dramatic difference in f opt I f x-ray between LIS and FRED outbursts ) as well as type I X-ray 
bursts suggests that a significant amount of mass is being accreted by the outer accretion 
disk and fed to the inner region. 



Radio detection of the source during LIS (JRupen et al.ll2004l ) with flat spectral slope 
hints toward a jet outflow. However, since the evolution of the source on OIR color-magnitude 
diagram is consistent with irradiation, it can be concluded that the jet, if present, is either 
quite weak in terms of the input power, or the transition from optically thick to optically thin 
jet emission occurs at a wavelength much longer than J-band. Current models for jets from 
compact systems suggest that this might imply a large jet launching area or the jet base, or 
much weaker particle acceleration (JMarkoff et al.ll200ll ). Thus our observations support the 
current hypothesis that jet-luminosity is not s ignificant in most neutron star X-ray binary 
systems ( jMigliari fc Fenderl l2006l : iFenderl |2006| ) . 



A possible scenario for the LIS might he that an M that would ordinarily he high enough 
to move the disk inwards, and thus trigger a soft state, somehow doesn't do that, and launches 
an outflow instead. Thus the standard processes associated with the DIM are disrupted by 
the absence of an inner disk in this high M state. If the accretion fiow in the inner region 
in a LIS is radiatively ine fficient, it can drive the obseryed outflow , as suggested b y recent 



theoretical works (see e.g. iBlandford fc Begelmarull999l : lMeierll200ll : lNarayanll2005l ). 
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Table 1. Date of peak flux during the outbursts of Aquila X-1 between MJD 50984 - 

MJD 54411 (1998 June 20 - 2007 Nov 7). 
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Fig. 1. — Multi- wavelength long-term monitoring of Aquila X-1 between MJD 50984 - MJD 
54411 (1998 June 20 - 2007 Nov 7). Top panel: Near-IR J-band light curve. Middle panel: 
Optical R-band light curve. The numbers in parentheses refer to the serial number assigned 
to the major outbursts in Table [H Bottom panel: Logarithm of soft X-ray (~1.5-12 keV) 
count-rate as measured by the All-Sky Monitor onboard the RXTE. Since the ASM does 
not detect the source when in quiescence, only those ASM data are plotted for which the 
signal-to-noise ratio is > 4. The R and J-band fluxes include the flux from Aquila X-1 as 
well as from the contaminating star. During periods of quiescence, the soft X-ray flux is 
below the detection limit of ASM but the combined flux from Aquila X-H- contaminating 
star remains detectable in both R and J bands. Note that a major outburst which occurred 
during late 2005 - early 2006 could be seen only with the ASM but not the ground based 
telescopes due to close angular proximity to the Sun. Therefore we have not labelled this 
outburst in the middle panel. 
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Fig. 2. — Zooming into the R-band and ASM light curves shown in Fig. [T], showing the 
detailed light curve morphology during major outbursts 1-3. FRED outbursts are marked 
by filled circles, LIS by open circles. In this figure we are using a linear scale for the ASM 
count-rates to emphasize the difference between LIS and FRED. 
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Fig. 3. — Same as Fig. [2], but for outbursts 4-6. 



-21- 



16.0 

oi 17.0 

E, 

OC 18.0 

19.0 


; • 
■ . 


• 
1 


• • •• . 

1 1 , , , , 


G 


Outburst #7 ' 

Q - 
O 

®0 G ^ 
. . 1 ' 


8.0 

~ 6.0 

o 

S 4.0 

CO 

^ 2.0 




>»« 


• 
• 


(fCD $ 




* : 
■ . 1 ■ 



53460 53470 



CD 

E, 


18.0 


DC 






19.0 




5.0 


'ST 

Q. 


4.0 
3.0 


CO 
< 


2.0 



1.0 



53940 



53480 



53490 
MJD 



53500 53510 



53520 



- 


• 


• 

• 
.... 1 


• 

• • 

1 1 


ffl ffl 
1 


Outburst #8 

G 
G O 




- 






.... 1 












53945 




53970 



54220 54230 54240 



54250 54260 
MJD 



54270 54280 54290 



Fig. 4. — Same as Fig. [2], but for outbursts 7-9. 
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Fig. 5. — Same as Fig. [2], but for outburst 10. 
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Fig. 6. — Optical-soft-X-ray correlation in Aquila X-1. Data from all the major out- 
bursts from MJD 50984 - MJD 54411 (1998 June 20 - 2007 Nov 7) for which the X-ray 
state c o uld be determined from contemporaneous RXTE observations ( Maccarone fc Coppi 



Yu fc Dolence 2007 



2003al lbl: iReig et al.l[20oi : iMaitra fc Bailvnll2004J : [Rodriguez et~aDl2006 
also Manuel Linares, priv. comm.). The dashed lines LI and L2 which follow the equations 
fa = 1.25fASM — 75 and //? = 1.5/a5a/ + 270 respectively, roughly form an envelope con- 
taining most of the FRED points. The dotted lines L3 and L4 which follow the equations 
fn = 9/yi5A/ — 240 and f^ = 12/^15^ + 240 respectively, roughly form an envelope containing 
most of the LIS points. Filled circles represent data taken during soft X-ray state and open 
squares are those during hard X-ray state. 
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Fig. 7. — Evolution of J-R color during outbursts 2,3,5 and 6 of Aquila X-1. The increase 
in brightness from quiescence is shown in the top panels, with the R-band data shown by 
filled circles and J-band data shown by open squares. The outburst serial number (Table. [1]) 
is given in parentheses. Variation of corresponding J-R color is shown in the bottom panels. 
Epochs of LIS as inferred from Figs. [THH] are labelled. Unless labelled as a LIS, the OIR/soft- 
X-ray light curve morphology resembles that of FRED outbursts. 
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Fig. 8. — Same as Fig. [71 but for outbursts 7-10. 
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Fig. 9. — Evolution of outbursts 2, 3, 5 and 6 of Aquila X-1 on the OIR color-magnitude 
diagram. Observations when the source was in hard/power law state during the rise of an 
outburst are shown by (A). The decaying hard state, are shown by (v)- Outbursts where 
state transitions were not reported are shown by (D). Observations during soft/thermal 
dominated state are shown by (o). The quiescent level is shown by the dashed line. A 
typical error bar is shown in the top left panel near (1.4, 17.75). The solid, curve from 
bottom-right to top-left in each panel shows the color and brightness expected by heating 
a disk-shaped blackbody of radius 5 light seconds, inclined to the observer at 45° and at a 
distance of 3 kpc. The quantity / represents a multiplicative area correction factor. The 
expected color and magnitude of the disk between temperatures of 4000 -6000 K are labelled. 
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Fig. 10. — Same as Fig. [9l but for outbursts 7, 8, 9 and 10. 



